Magnetic phase diagram of MnSi inferred from magnetization and ac susceptibility 
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We report simultaneous measurements of the magnetization and the ac susceptibility across the 
magnetic phase diagram of single-crystal MnSi. In our study we explore the importance of the 
excitation frequency, excitation amplitude, sample shape, and crystallographic orientation. The 
susceptibility, fj,o dAf/dB, calculated from the magnetization, is dominated by pronounced maxima 
at the transition from the helical to the conical and the conical to the skyrmion lattice phase. The 
maxima in dM/dB are not tracked by the ac susceptibility, which in addition varies sensitively 
with the excitation amplitude and frequency at the transition from the conical to the skyrmion lattice 
phase. The same differences between /io dM/d£> and the ac susceptibility exist for Mni-^Fe^Si (x = 
0.04) and Fei-^Co^Si [x = 0.20). Taken together our study establishes consistently for all major 
crystallographic directions the existence of a single pocket of the skyrmion lattice phase in MnSi, 
suggestive of a universal characteristic of all B20 transition metal compounds with helimagnetic 
order. 

PACS numbers: 75.30.Kz, 75.30.Gw, 75.25.-j 



I. INTRODUCTION 

B20 transition metal compounds such as MnSi, 
Fei-zCo^Si, or FeGe exhibit long-wavelength helimag- 
netic order at low temperatures. As emphasized already 
in Landau-Lifshitz, Vol. VIII, Sec. 52jHthe helimagnetic 
order in B20 compounds reflects a well defined set of hier- 
archical energy scales. On the strongest scale ferromag- 
netic exchange interactions favor parallel spin- alignment. 
Due to the lack of inversion symmetry of the B20 crys- 
tal structure this is followed on intermediate scales by 
Dzyaloshinsky-Moriya (DM) spin-orbit interactions, fa- 
voring perpendicular spin alignment. The combination 
of ferromagnetic exchange and DM interactions gener- 
ates a long wavelength helical modulation. The chirality 
of the DM interaction and thus of the helical modulation 
is thereby determined by the enantiomorph of the B20 
structure. Finally, on the weakest energy scale higher or- 
der spin-orbit coupling, also referred to as crystal electric 
field effects, aligns the propagation direction of the he- 
lical modulations along either (111) or (100) depending 
on the system. 

Recently, interest in the spin order in B20 transition 
metal compounds has revived, when a small phase pocket 
just below the helimagnetic transition temperature T c , 
long known as the A phase, was identified as the first 
example of a skyrmion lattice in a magnetic materialPEl 
The skyrmion lattice represents a novel form of magnetic 
order composed of spin vortices with a non-zero topo- 
logical winding number. The skyrmion lattice is thereby 
always perpendicular to the applied field; i.e., the pinning 
with respect to the crystal lattice is extremely weakPEl 
Perhaps most remarkable, the coupling of electric cur- 
rents to the skyrmion lattice is extremely efficient, lead- 
ing to emergent electrodynamics an d siz able spin torque 
effects at ultra-low current densities!^ 

The microscopic identification of the skyrmion lat- 
tice was initially based on small-angle neutron scatter- 



ing (SANS) in MnSi, Mni_ x Fe x Si, Mni^^Co^Si, and 
Fei-zCo^SiPE] The non-zero topological winding num- 
ber, as the defining aspect of the skyrmion lattice, was 
first inferred from topological contributions to the Hall 
effeclP and recently confirmed in a rather sophisticated 
small-angle neutron scattering study™ Real-space imag- 
ing by means of Lorentz force microscopy confirmed 
the formation of the skyrmion lattice in thin samples 
of selected B20 transition metal compounds, notably 
Fei-zCo^Si, FeGe, and MnSiPEU As the most recent 
development a skyrmion lattice has been observed in 
the multiferroic compound C^OSeOaj 13 * 14 ! which crys- 
tallizes in the same P2i3 space group as the B20 com- 
pounds. 

The large body of macroscopic and microscopic infor- 
mation on the magnetic phase diagram and skyrmion 
lattice phase in B20 transition metal bulk compounds 
available so far has been explained in terms of Ginzburg- 
Landau theory taking into account thermal Gaussian 
fluctuations.^ For the thin samples studied by Lorentz 
force microscopy the effects of reduced spatial extent were 
considered, as they destabilize the conical phase, thereby 
supporting other routes to skyrmion lattices.^ Neither 
the data of microscopic properties nor these theoretical 
models suggest the existence of a more complex phase 
diagram with multiple phases. 

In contrast, thirty years ago Kadowaki et al. reported 
that the A phase in MnSi consists of two phase pockets,^ 
putatively suggesting several phases with different spin 
structures of nontrivial topology. The evidence for the 
two phase pockets was thereby purely based on the pres- 
ence of two kinks in the temperature dependence of the 
magnetization when traversing the A phase near its up- 
per and lower boundary. In contrast, no kinks were ob- 
served in the magnetization for a field in the center of 
the A phase, suggesting the absence of an A phase. How- 
ever, Kadowaki et al. recorded the magnetization with 
a Faraday balance which requires a field gradient. In 
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turn, the magnetization data shown by Kadowaki et al. 
may have been subject to parasitic magnetic torques and 
did not represent the longitudinal magnetization. More- 
over, as pointed out by Kadowaki et al. the temperature 
control of their setup was rather unsatisfactory. Data 
were therefore recorded by means of temperature sweeps. 
This may have caused small systematic temperature dif- 
ferences between different temperature sweeps. Finally, 
magnetoresistance data as a function of magnetic field, 
also reported by Kadowaki et al., suggested the presence 
of only a single pocket of the A phase. 

The question of several phase pockets was recently fol- 
lowed up in a study of FeGe, where an even more complex 
magnetic phase diagram including claims of mesophases 
with liquid-crystal-like appearance were proposed.^ Yet, 
the phase diagram in FeGe was purely inferred from the 
ac susceptibility as measured at a unique excitation fre- 
quency of 1 kHz in a Quantum Design PPMS. Moreover, 
the FeGe sample studied, being vapor transport grown, 
had an irregular shape causing ill-defined demagnetizing 
fields and a poor residual resistivity ratio characteristic 
for defects and disorder (see also the note added at the 
end of this paper concerning C^OSeOa). 

The possible existence of multiple pockets in the 
magnetic phase diagram of MnSi and FeGe has been 
explained in terms of a Landau theory, where the 
ground state is stabilized by a combination of magnetic 
anisotropics and a controv ersial gradient term favoring 
soft amplitude fluctuations! 17 * 18 ! The possible existence 
of a more complex magnetic phase diagram is therefore 
believed to provide putative evidence for a different mi- 
croscopic mechanism at the heart of the skyrmion lattice 
than that considered so far. 

To clarify the magnetic phase diagram in MnSi and 
other B20 transition metal compounds we report in the 
following simultaneous measurements of the magnetiza- 
tion and ac susceptibility in high-quality single-crystal 
MnSi. In our study we carefully explored the impor- 
tance of the excitation frequency and amplitude for the 
ac susceptibility, as well as the sample shape and crys- 
tallographic orientation. As the skyrmion lattice phase 
is embedded in a spin-flop phase, also known as conical 
phase, we covered the entire magnetic phase diagram. 
We finally also performed exploratory measurements in 
Mni_ x Fe x Si and Fei-^Co^Si to demonstrate the generic 
nature of our observations for other B20 transition metal 
compounds and in the presence of disorder and defects. 

In our study of MnSi we establish consistently for all 
major crystallographic directions the existence of a sin- 
gle pocket of the A phase. In contrast, there is no evi- 
dence whatsoever pointing at any additional subphases. 
The most prominent aspect of our results concerns the 
observation of maxima in the susceptibility fiQ dM/dB, 
calculated from the magnetization, at the border be- 
tween the helical, the conical and the skyrmion lattice 
phase. General thermodynamic considerations and the 
imaginary part of the ac susceptibility establish, that 
the maxima in fiQ dM/dB are the characteristic of pro- 



cesses on extremely slow time-scales involving very large 
magnetic objects, consistent with a crossover or slightly 
broadened first-order transitions at the helical to conical 
transition or the boundary of the A phase, respectively. 
Surprisingly, the maxima in (iq dM/dB are not tracked 
by the real part of the ac susceptibility, Rexac, which 
even varies sensitively with the excitation amplitude and 
frequency at the transition to the skyrmion lattice. Re- 
producing the temperature dependence of the magnetiza- 
tion reported by Kadowaki et al. (Ref.^), we are finally 
able to show that the magnetic phase diagram reported 
in Ref.QSlis not correct. 

Moreover, our results question the evidence for the 
rather complex magnetic phase diagram of FeGe reported 
by Wilhclm et al.pl which was purely based on the ac 
susceptibility recorded at a single frequency of 1 kHz. 
Taken together with data we recorded for Mn^^FezSi 
and Fei-^CozSi, our study suggests that the magnetic 
phase diagram we observe in MnSi represents a univer- 
sal characteristic of all B20 transition metal compounds 
with helimagnetic order. 



II. EXPERIMENTAL METHODS 

For our studies of the magnetic phase diagram of 
MnSi a large single crystal was grown in a bespoke ul- 
trahigh vacuum compatible image furnace!^ The feed 
rods and seed rods used for the float zoning were pre- 
pared from high-purity elements and cast in a UHV- 
compatiblc radio-frequency induction furnace equipped 
with a bespoke all-metal-sealed Huykin crucible and cast- 
ing moldPS Samples were oriented by means of Laue x- 
ray diffraction and cut with a wire saw. The sample 
surfaces were gently polished. The residual resistivity 
ratio (RRR) of samples prepared from the same ingot 
was very high exceeding 300, characteristic of an ex- 
cellent sample quality. For the preliminary studies in 
Mn 1 _ 2 .Fe :E Si and Fei-^Co^Si reported at the end of this 
paper single-crystal samples were grown following the 
same procedure as for MnSi. Previous studies on the 
same Mni_a;Fe x Si an d Fei-^Co^Si single crystals were 
reported in Refs.^il. 

To track the importance of demagnetizing fields and 
crystallographic orientation on the magnetic phase dia- 
gram of MnSi four samples with three different sample 
geometries were studied as summarized in Table[I] The 
sample shapes investigated consisted of a bar, two cubes, 
and a thin platelet. Most of our measurements were car- 
ried out on the bar-shaped single-crystal denoted as sam- 
ple 1. The dimensions of sample 1 were 6 x 1 x 1mm 3 
with the long direction parallel to a (100) axis. The small 
faces of sample 1 were perpendicular to (110) axes. The 
magnetic field was applied along the long direction of 
sample 1 to minimize the effects of demagnetizing fields. 
The samples of Mni-^Fe^Si and Fei-^Co^Si, denoted 
as samples 5 and 6 in Table[TJ had the same shape and 
orientation as sample 1. 
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TABLE I: Samples investigated in this study. Magnetic fields were applied along the direction labeled (a). 



system sample shape (a) size, orientation (b) size, orientation (c) size, orientation N 



MnSi 


1 


bar 6 mm, (100) 


1 mm, 


(110) 


1 mm, 


(110) 


0.074 


MnSi 


2 


cube 1 mm, (100) 


1 mm, 


(110) 


1 mm, 


(110) 


0.333 


MnSi 


3 


cube 1 mm, (110) 


1 mm, 


(111) 


1 mm, 


(211) 


0.333 


MnSi 


4 


platelet 0.115 mm, (100) 


3 mm, 


(110) 


1 mm, 


(110) 


0.837 


Mni-aFe^Si, 


x = 0.04 5 


bar 6 mm, (100) 


1 mm, 


(110) 


1 mm, 


(110) 


0.074 


Fei-^CozSi, 


a; = 0.20 6 


bar 6 mm, (100) 


1 mm, 


(110) 


1 mm, 


(110) 


0.074 



Samples 2 and 3 were cubic and served to study 
the role of the crystallographic orientation under un- 
changed sample geometry. Sample 2 was oriented such 
that two opposing faces were oriented perpendicular to 
a (100) axis. The other four faces were oriented per- 
pendicular to a (110) axes. This contrasts sample 3, 
which was oriented such that the faces were perpendicu- 
lar to either (111), (110), or (211), respectively. Sam- 
ples 2 and 3 allowed us to reproduce data for (110) 
and by rotating the cube through 90° to record data 
for field parallel to a (100) or (111) under unchanged 
sample geometry. Sample 4 was a thin platelet of 
3 x 1 x 0.115mm 3 (length x width x thickness); i.e., the 
dimensions were characteristic of samples used in mag- 
netotransport measurements. The smallest direction of 
the platelet was thereby parallel to a (100) axis, while 
the long direction, equivalent to the current direction in 
transport measurements, was parallel to a (110) axis. 

The demagnetizing factors of the samples were deter- 
mined by approximating the sample shape as a rectan- 
gular prism. For the field directions given above they 
are Ni = 0.074, N 2 = N 3 = 0.333, and N 4 = 0.837 for 
samples 1, 2, 3, and 4, respectively. As summarized be- 
low all temperature and field dependencies were found 
to be qualitatively independent of the sample geometry. 
Demagnetizing fields only affected the absolute values of 
the transition fields in excellent agreement with the cal- 
culated demagnetization factors given above. In addi- 
tion, the quantitative size of the anomalies at the phase 
boundaries were found to decrease with increasing de- 
magnetizing factor while the range of the inconsistencies 
increased. For lack of space, all experimental data are 
shown as a function of applied field without correction for 
demagnetizing fields, while the phase diagrams inferred 
from the data are shown as a function of both applied 
and internal field, respectively. 

The magnetization and the ac susceptibility were mea- 
sured in a Quantum Design physical properties measure- 
ment system (PPMS). Data points were recorded at fixed 
temperature and magnetic field values. Unlike common 
practice in many studies no data were recorded while 
the field or the temperature were swept continuously. At 
each given temperature and field value the magnetiza- 
tion was determined by means of a standard extraction 
method. This was followed by the measurement of the 



ac susceptibility. The ac field was only turned on while 
recording the susceptibility data. We note that all data 
reported here were recorded using the same sample ther- 
mometer. However, in comparison to other studies re- 
ported in the literature there may be small systematic 
differences of the absolute temperature values on the or- 
der of a few % (see also the clarification of Ref.^). The ac 
susceptibility was recorded in the frequency range from 
10 Hz to 10 kHz. For typical values of the resistivity and 
the susceptibility in the regime of the A phase of MnSi, 
PMnSi(28K) w 35 /fJJcm, XMnSi ~ 0-3, and the skin depth 
<5 M „Si(28K) = (2p/(c^o(x + l))) 1/2 « 0.26-/- 1 /2 ms -i/2. 
In the range of excitation frequencies from 10 Hz to 
10 kHz applied in our study the skin depth was between 
82 mm and 2.6 mm, respectively, and thus considerably 
larger than the relevant sample dimensions. As the resis- 
tivity of Mni^Fe^Si and Fei^CozSi is larger than that 
of MnSi the same conclusion is true for these compounds. 

In the course of our studies we became aware of two 
minor technical limitations of the PPMS as concerns si- 
multaneous measurements of the ac susceptibility and 
magnetization. First, when repeating measurements of 
the ac susceptibility on the same sample and sample sup- 
port (remounting the sample from scratch between mea- 
surements), the signal may include a small temperature- 
independent signal contribution that varies slightly be- 
tween different measurements. Second, the ac suscepti- 
bility measured for parameters where highly reversible 
behavior is expected, i.e., at high fields and/or high tem- 
peratures, differs by several % from the susceptibility, 
HodM/dB, calculated from the magnetic field depen- 
dence of the magnetization. 

We have therefore recalibrated the measured ac sus- 
ceptibility of all samples and sample orientations with 
respect to the susceptibility hq&M/AB calculated from 
the magnetization, by the same procedure as follows. We 
subtracted the same constant value of all ac susceptibil- 
ity data, Ax = 0.025, at all temperatures and fields. 
We then multiplied all of our ac susceptibility data by 
a factor of 1.1. This way the ac susceptibility and val- 
ues of no dM/dB agreed for field and temperature values 
well away from the magnetically ordered state, notably at 
high magnetic fields and/or high temperatures. We be- 
lieve that the small deviation visible in Fig.[3](a) reflects 
this mismatch. Taken together, the conclusions of the 
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work reported here are independent of this recalibration 
of the ac susceptibility data. 



III. EXPERIMENTAL RESULTS 

The presentation of our experimental results is orga- 
nized as follows. We begin with the magnetic field depen- 
dence of sample 1, where the ac susceptibility data were 
recorded at a low excitation amplitude of 0.1 mT and 
an excitation frequency of 10 Hz unless explicitly stated 
differently. This is followed by a comparison of typical 
changes of the ac susceptibility for larger excitation am- 
plitudes, notably 0.3 mT and 1 mT at 10 Hz, and higher 
excitation frequencies of 91 Hz and 911 Hz at 0.1 mT. 

In the second part of this section we present the tem- 
perature dependence of the magnetization and ac suscep- 
tibility of sample 1 for 0.1 mT and 10 Hz, followed by data 
illustrating the role of changes of excitation amplitudes 
and frequencies. This way we establish the consistency 
of the rather unusual temperature dependence with the 
magnetic field dependence. In particular, we are able to 
show the consistency of the evidence for a single A phase 
seen in the magnetic field dependence with the data re- 
ported by Kadowaki et aZ.P^ 

The third part of the experimental results section con- 
cerns variations of the ac susceptibility under changes of 
sample shape and crystallographic orientation. On the 
one hand this allows us to point out potential pitfalls 
when inferring the magnetic phase diagram of B20 com- 
pounds from ac susceptibility data alone. On the other 
hand our studies reveal a regime of inconsistency be- 
tween magnetization and ac susceptibility at the bound- 
ary of the A phase. The results section concludes with 
a presentation of preliminary results on Mni_ x Fe x Si and 
Fei-zCo^Si. 

Given the amount and complexity of the data pre- 
sented in the following it is helpful to present at first 
in Fig.[T] a qualitative summary as a guide of the essen- 
tial features of the magnetization in the regime of the 
skyrmion lattice phase. The key characteristic of the 
magnetization in the skyrmion lattice phase is a plateau, 
that is slightly tilted with respect to the surrounding 
phases. At the boundary of this plateau the magnetiza- 
tion changes discontinuously (first order). This discon- 
tinuous step is smeared out in the experimental data. As 
a function of increasing magnetic field, depicted by trace 
(i), the magnetization shows always two first-order tran- 
sitions at which the magnetization increases. However, 
as a function of increasing temperature more complex be- 
havior is observed. For a temperature scan in the vicin- 
ity of the upper boundary of the skyrmion lattice phase, 
depicted by trace (ii), the magnetization displays a de- 
creasing step followed by an increasing step. In contrast, 
in a temperature scan in the vicinity of the lower bound- 
ary, depicted by trace (iii), the magnetization shows an 
increasing step followed by a decreasing step. Finally, for 
an intermediate field it is, in principle, possible to obtain 




FIG. 1: (Color online) Schematic depiction of the magneti- 
zation as a function of temperature and magnetic field. In 
the regime of the skyrmion lattice phase the magnetization 
is characterized by a tilted plateaux. The boundary of this 
plateaux is a discontinuous (first order) step. The traces il- 
lustrate typical scans as a function of magnetic field or tem- 
perature. In the experimental data the discontinuous steps 
are smeared out. 



temperature scans that show zero or one transitions, cf. 
trace (iv). This is a purely accidental effect, which ac- 
counts for our data as well as that reported by Kadowaki 
et aZ.P^Yet, behavior such as that depicted by trace (iv) 
is clearly no evidence of the existence of several phase 
pockets. 

As part of the presentation of our experimental results 
we define various characteristic transition and crossover 
fields and temperatures. All definitions arc thereby based 
on /j, dM/dB and [1qAM/AB with corresponding fea- 
tures in Rex ac and Imxac- Characteristic field values 
were determined under increasing magnetic fields. Unfor- 
tunately even for a minimalist description it is necessary 
to account for four transitions as a function of field (heli- 
cal — » conical — > skyrmion lattice — > conical — » ferromag- 
netic) and three transitions as a function of temperature 
(helical — > paramagnetic; conical —¥ skyrmion lattice — ¥ 
paramagnetic). When trying to track the width of the 
transitions this requires a large number of characteristic 
fields and temperatures. Even though it might appears 
tedious, we feel that this is necessary to demonstrate the 
full consistency between field and temperature depen- 
dent data. We thereby note, that some of the definitions 
given are rather empirical and cannot be followed in the 
data very accurately. For clarity we use a short-hand 
notation to mark all characteristic field and temperature 
values in the figures. A summary of how the field and 
temperature values are defined is shown in Tables |ll| and 
|III| More detailed definitions are given in the text, when 
the experimental evidence for the various transitions is 
presented first. As mentioned above all data are pre- 
sented as measured experimentally, i.e., as a function of 
applied magnetic field. The magnetic phase diagrams are 
shown as a function of applied magnetic field as well as 
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a function of internal magnetic field, where the effects of 
demagnetizing fields have been corrected. 



A. Magnetic field dependence 

Shown in Fig.[2](a) is the magnetization M of sample 
1 as a function of magnetic field along (fOO) for temper- 
atures in the vicinity of the helimagnetic transition tem- 
perature T c . We begin with an overview and an account 
of the global features and return to a detailed comparison 
including the definition of the transition fields in Fig. [3] 
Data for both increasing and decreasing fields are shown, 
where the field sweeps were started at — 1 T and +1 T for 
increasing and decreasing fields, respectively. For this 
field history data correspond to the field-cooled behavior 
in temperature scans presented below. None of the field 
sweeps displayed hysteresis between increasing and de- 
creasing field at the resolution of the methods used. All 
our data are perfectly consistent with previous studies. 

Starting from B = the magnetization increases al- 
most linearly. Well below T c the initial quasi-linear in- 
crease for B < B c i arises from an increasing anharmonic- 
ity of the helical modulation, a slight reorientation of the 
helical modulation, and changes of domain population as 
seen in neutron scattering studies; see, e.g., Refs. 23 ~ 25 . 
Depending on the field direction the reorientation at B c \ 
represents a cross-over or a symmetry breaking second 
order phase transition. In contrast, for B > B c \ the 
linear increase is a characteristic of a closing-in of the 
spin- flop phase (the conical phase) up to B C 2, the tran- 
sition to the field-polarized ferromagnetic state. On the 
field scale shown in Fig. [2] the magnetization appears al- 
most saturated above B C 2- However, as reported in the 
literature up to ~33 Tf^Sl the highest fields studied, the 
magnetization is unsaturated and nonlinear characteris- 
tic of itinerant magnetism. On the field scale shown here, 
the rather small changes of the magnetization in the A 
phase are not obvious. 

Deviations from the quasilinear field dependence of the 
magnetization are best seen in the susceptibility. The 
real part of the ac susceptibility Re^ac, is shown in 
Fig.|2j(b) for the same field and temperature values as 
the magnetization shown in Fig.[2](a). At low tempera- 
tures and low fields the ac susceptibility increases gently 
before reaching a well-defined plateau above B c \ . Again 
data for increasing and decreasing fields are shown, where 
we find no hysteresis between increasing and decreasing 
fields. 

For temperatures just below T c a reduced value of 
the susceptibility is observed for magnetic fields between 
~ 100 mT and ~ 200 mT. Previous studi es hav e estab- 
lished that this corresponds to the A phasePZESln other 
words the magnetization in the A phase is characterized 
by a plateau of reduced slope as a function of increas- 
ing field, which implies that there must be a regime of 
increased slope connecting the conical phase with the A 
phase at its boundaries. Indeed, as shown in Fig.[2](b) 
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FIG. 2: (Color online) Typical magnetization and ac suscepti- 
bility data of sample 1 as a function of applied magnetic field. 
Data were recorded for various temperatures in the regime 
of the A phase. The ac susceptibility data were recorded for 
an excitation amplitude of O.f mT parallel to (100) at 10 Hz. 
See also Fig.[3]for definitions of the transition fields, (a) Mag- 
netization as a function of applied magnetic field, (b) Real 
part of the ac susceptibility as a function of applied magnetic 
field, (c) Imaginary part of the ac susceptibility as a function 
of applied magnetic field. Data shown in panel (c) have been 
shifted by constant value of 0.025 for clarity. 



for an excitation amplitude of 0.1 mT and excitation fre- 
quency of 10 Hz chosen here, Re Xac displays narrow max- 
ima at the transition between the conical phase and the 
A phase. 

The imaginary part of the ac susceptibility, Imx ac , 
shown in Fig.[2](c), shows clear maxima at the transition 
fields between the conical phase and the A phase [note 
that data shown in Fig.[2](c) have been shifted by a con- 
stant value of 0.025 for clarity]. This provides evidence 
that the narrow maxima in Re \&c at the border of the A 
phase reflects the response of large magnetic objects on 
very slow time scales which originate in a slightly broad- 
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TABLE II: Summary of characteristic transition and cross-over fields and associated labels used in the figures shown in this 
paper. The signatures in Ho dM/dB were used to define the actual field values; these correspond within the accuracy of our 
data to the signatures in ReXac and ImXac listed in this table. The definition of the signatures refers to data recorded for 
increasing magnetic field. Note that our definition of B^i an d B^ 2 corresponds to values referred to as Bai and Ba2 in the 
literature. 



in text in figures signature in Ho dM/dB 

B~ (1) begin /LtodM/dB^Rexac 

Bci (2) peak 

B+ (3) end ModM/dB/ReXac 

(4) begin of peak 

B^ (5) begin of plateau of reduced 

Ho dM/dB in A phase 

^A2 (6) end of plateau of reduced 

Ho dM/dB in A phase 

B^ 2 (7) end of peak 

B C 2 (8) point of inflection 



signature in Im x a 



signature in Re Xac 

begin Ho dM/dB^Re x? 
point of inflection 
end /^odM/dB/ReXac 
begin of deviation from 
Re Xac in conical phase 
begin of plateau of reduced end peak 
Re Xac in A phase 
end of plateau of reduced 
Re Xac in A phase 
end of deviation from 
Re Xac in conical phase 
point of inflection no signature 



no signature 
peak for large amplitude 
no signature 
begin peak 



begin peak 
end peak 



TABLE III: Summary of characteristic transition and cross-over temperatures and associated labels used in the figures shown 
in this paper. The signatures in hoAM/AB were used to define the actual temperature values; these correspond within the 
accuracy of our data to the signatures in ReXac and ImXac listed in this table. The definition of the signatures refers to data 
recorded for increasing temperature. Note that our definition of T^~i an d T£ 2 corresponds to values referred to as Tai and Ta2 
in the literature. 

in text in figures signature in Ho^-M / AB signature in ReXac signature in ImXac 

T^j (i) begin of maximum begin of deviation from begin peak 

Re Xac in conical phase 
(ii) end of maximum end /xoAAf/AB^ReXac end peak 

T^ 2 (iii) end of deviation from end of deviation from no signature 

HqAM/AB in conical phase ReXac in conical phase 
T C 2 (iv) point of inflection point of inflection no signature 



ened first-order step in the magnetization, i.e., broad- 
ened S anomalies in ReXac- The finite imaginary part 
of the ac susceptibility provides unambiguous evidence 
of the dissipation typically observed in regimes of very 
slow dissipative processes. This does not contrast the 
absence of hysteresis between increasing and decreasing 
fields within the resolution of the methods used, since the 
pinning may just be very weak. 

Neither Re Xac nor Im Xac display any additional struc- 
ture such as double peaks or additional shoulders at the 
phase boundary or deep inside the A phase. Thus there 
is no evidence that would hint at additional phase transi- 
tions. This is consistent with a wide range of microscopic 
studies we have performed (mostly neutron scattering), 
which do not provide any experimental evidence suggest- 
ing several pockets of the A phase either. 

It is now instructive to compare in further detail 
the susceptibility calculated from the magnetization, 
Ho dM/dB, shown in Fig.[2](a) with the real part of the ac 



susceptibility, ReXac shown in Fig.[2](b). At T = 28.5 K 
shown in Fig.[3](a), i.e., above the helimagnetic transi- 
tion, Ho dM/dB and Re x ac agree essentially (we believe 
that the small difference of the data is due to a small 
systematic error as explained above). 

With decreasing temperature, shown in Figs.[3](b) 
through |3](f), we find in the regime of the helical to 
conical transition that Ho dM / dB increases much faster 
than Re Xac and displays a pronounced maximum and an 
abrupt drop. The absence of hysteresis between increas- 
ing and decreasing fields is thereby characteristic of an 
extremely soft magnetic state in which the ac suscepti- 
bility, probing tiny hysteresis loops, does not fully track 
Ho dM/dB. At the same time we find a tiny contribu- 
tion to ImXac within our detection limit. This suggests 
that hysteretic losses are very small at the helimagnetic 
to conical transition and reveals a peculiar characteristic 
of the magnetic properties of MnSi, which, to the best of 
our knowledge, has not been noticed before. 
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FIG. 3: (Color online) Comparison of the magnetic field dependence of the magnetization, M, the susceptibility calculated 
from the magnetization, /io dM/dB, and the real part of the ac susceptibility, Rexac- Data shown here were recorded in sample 
1 for various temperatures in the regime of the A phase. The ac susceptibility data were recorded for an excitation amplitude 
of 0.1 mT along (100) at 10 Hz. 



To be able to track the characteristics of the helical 
to conical transition wc define three field values (see also 
Tablepl|. The lower field, B~ 1: denotes the point where 
/io dM/dB and Re Xac begin to differ, the transition field 
B c i is located at the maximum of /io dM/dB, and finally 
there is the upper field B^ 1: above which /io dM/dB and 
R- e Xac no longer differ. Values of £? cl agree thereby with 
the transition fields reported in the literature. The sim- 
plest explanation for the deviation between /io dM/dB 
and Re Xac is the response of large magnetic objects at 
very slow time scales for instance due to a reorientation 
of the helical modulation. 

We turn now to the A phase, where the same quali- 
tative difference between /i dM/dB and Rexac emerges 
like that seen near B c \\ see also Figs.|3](b) through [3](f). 
Here, clear maxima exist in /io dM/dB. Yet, these max- 
ima are not tracked by Rexac which is always smaller 
than /i dM/dB. This similarity between the boundary 
of the A phase and the helical to conical transition, which 
arise microscopically on different ground, implies that the 
phenomenology in the magnetization and ac susceptibil- 
ity by itself cannot be taken as evidence for the formation 
of mesophases of nontrivial spin topology. 

Analogous to the helical to conical transition we de- 
fine the following transition fields of the A phase. B^ 
denotes the onset of the maximum in /to dM/dB at the 
transition to the A phase. For most of the parame- 
ter range studied this agrees with the field value where 
/i dM/dB and Rexac differ first. Likewise B~^ x repre- 
sents the upper field where the maximum in /io dM/dB 
is completed; in the ac susceptibility this corresponds to 
the point where /io dM / dB and Re Xac merge and the 
plateau of the reduced ac susceptibility in the A phase 
begins. Analogous transition fields B~^ 2 and B\ 2 are de- 



fined at the second maximum in /io dM/dB and the cor- 
responding features in Rexac, respectively. By compari- 
son with previous work BJ^ X corresponds to the low-field 
boundary reported in the literature, typically denoted as 
£?ai, while B^ 2 corresponds to the high-field boundary 
of the A phase, typically denoted as Ba2- 

A key observation of our study concerns the absence of 
the narrow maxima in Rexac for sufficiently large exci- 
tation frequencies as compared with /io dM/dB for most 
of the temperatures studied. In fact, Rexac frequently 
displays a reduction, where the onset of the reduction 
agrees with B^ x in /i dM/dB. The differences between 
/io dM/ dB and Re Xac thereby depends on the precise 
choice of the excitation amplitude and frequency as illus- 
trated in Figs. [4] and [5j 

The importance of the excitation frequency is illus- 
trated in Fig.|2](a), which displays data at two typical 
temperatures in the A phase for various frequencies and 
the same excitation amplitude of 0.1 mT. Here the nar- 
row maxima at the border of the A phase decrease with 
increasing excitation frequency. As emphasized above, 
this has nothing to do with the skin depth. For the stan- 
dard excitation frequency of 911 Hz, typically used by the 
PPMS, no maxima are observed. This compares with the 
amplitude dependence recorded for a fixed excitation fre- 
quency of 10 Hz as shown in Figs.[5](a) and[5](b). Here 
the narrow maxima increase with increasing amplitude. 

Even though the precise features in the ac suscepti- 
bility at the border of the A phase vary and depend on 
the measurement parameters, the values of B^ and B\ 2 
arc essentially unchanged; i.e., the outer boundary of the 
A phase in most studies reported in the literature were 
determined correctly. However, even though there is a 
pronounced frequency and amplitude dependence neither 
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FIG. 4: (Color online) Typical magnetic field depen- 
dence of the susceptibility inferred from the magnetization, 
fio dM/dB, and the real and imaginary parts of the ac suscep- 
tibility, Re Xac and Im Xac , respectively, at selected tempera- 
tures for various excitation frequencies. Data were recorded 
for an excitation amplitude of 0.1 mT. Curves at different tem- 
peratures have been shifted by constant values for clarity. See 
Table [IT] for a summary of the definitions of characteristic field 
values. 



FIG. 5: (Color online) Typical magnetic field depen- 
dence of the susceptibility inferred from the magnetization, 
fj,o dM/dB, and the real and imaginary parts of the ac suscep- 
tibility, Re Xac and Im^ac, respectively, at selected tempera- 
tures for various excitation amplitudes. Data were recorded 
for an excitation frequency of 10 Hz. Curves at different tem- 
peratures have been shifted by constant values for clarity. See 
Table |n] for a summary of the definitions of characteristic field 
values. 



the real nor the imaginary part of the ac susceptibility 
displays additional structure such as a double peak or an 
additional shoulder at the boundary or deep inside the 
A phase that might hint at additional phase transitions. 
Instead they are characteristic of a small range with a 
response of large magnetic objects on slow time scales, 
consistentwith a slightly smeared first-order transition. 
Thus, the frequency and amplitude dependence under- 
score an inherent risk of inferring erroneous features of 
the phase diagram when the ac susceptibility is not mea- 
sured carefully. 



B. Temperature dependence 

We now turn to the temperature dependence of the 
magnetization and the ac susceptibility to establish 
the consistency with the field-dependent data reported 
above. This brings us also to the origin of the phase dia- 
gram reported by Kadowaki et al.W^It is thereby helpful 
to keep in mind the qualitative picture shown in Fig.[T] 
where the discontinuous steps in the the magnetization 
are smeared out in the experimental data. 



Shown in Fig.[6](a) is an overview of the temperature 
dependence of the magnetization of sample 1 in the vicin- 
ity of the helimagnetic to paramagnetic transition for 
various applied magnetic fields up to 225 mT. The ac 
susceptibility data shown in Fig. [6] were recorded for an 
excitation amplitude of 0.1 mT at 10 Hz along (100). In 
order to determine the temperature dependence of the 
susceptibility as inferred from the measured magnetiza- 
tion, we have measured M(T) at field values that differed 
by 5mT and computed [i AM/ AB as a reasonable ap- 
proximation of fiQ dM/dB. 

We begin by describing the data qualitatively and re- 
turn to detailed definitions of characteristic transition 
and cross-over temperatures further below. Again a 
shorthand notation is used for the characteristic temper- 
atures, which is summarized in Tablc |Hl] 

With increasing temperature the magnetization ini- 
tially varies only weakly. For temperatures T > T c a 
point of inflection marks the onset of a strong Curie Weiss 
dependence as a key characteristic of the paramagnetic 
state at high temperatures. Details of the transition re- 
gion are best illustrated in the ac susceptibility as shown 
in Fig. [61(b). For instance, the point of inflection seen 
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FIG. 6: (Color online) Typical temperature dependence of the 
magnetization and the ac susceptibility of sample 1 at various 
magnetic fields in the regime of the A phase. Data shown here 
were recorded in sample 1. The ac susceptibility data were 
recorded for an excitation amplitude of 0.1 mT along (100) at 
10 Hz. See Fig.[7|for definitions of the transition temperatures, 
(a) Magnetization as a function of temperature, (b) Re^ac 
as a function of temperature, (c) Im^ac as a function of 
temperature. Data in panel (c) have been shifted by constant 
values of 0.025 for clarity. 



in the magnetization at a temperature slightly above the 
helimagnetic transition is also seen in Rex ac . In finite 
magnetic fields this point of inflection remains unchanged 
for fields up to ~ 100 mT and is shifted to lower temper- 
atures for slightly larger fields. As reported elsewhere, 
we have recently been able to show that the intermedi- 
ate regime between the helimagnetic transition and the 
point of inflection (marked IM in Figs. 12 and 20 1 exhibits 
all the characteristics of a Brazovskii transition 29 in con- 
trast to earlier claims of a skyrmion liquid or Bak- Jensen 
transitionPHU 

Between ~ 115 mT and ~ 210 mT, i.e., in the field 
range of the A phase, the most prominent feature of 



Rexac is a maximum along the upper and lower field 
boundary of the A phase that gives way to a minimum 
for intermediate fields. This is accompanied by a max- 
imum in Imx ac in a small temperature interval on the 
low-temperature side of the A phase as shown in Fig.[6](c) 
(for clarity curves have been shifted by a constant value 
of 0.025 in this panel). Marked by small arrows are in 
Fig.[6](c) the lower and upper boundary of the maximum 
in Imxac- As shown below this corresponds to the tran- 
sition temperatures and in [XqAM/ AB. 

A detailed comparison of the temperature dependence 
of no AM/ AB with Re%ac at selected magnetic fields is 
shown in Fig.[7j As in Fig. [6] the ac susceptibility data 
were recorded simultaneously for an excitation frequency 
of 10 Hz and an excitation amplitude of 0.1 mT along 
(100). In the following we define various transition and 
cross-over temperatures. As for the magnetic field de- 
pendence a main observation in our study concerns the 
differences between Rexac and hqAM/AB. 

Starting with the temperature dependence at 215 mT, 
shown in Fig.[7](a), the magnetization shows a shallow 
minimum that is accompanied by a pronounced maxi- 
mum in hqAM/AB in the same temperature range. This 
behavior corresponds to trace (ii) in Fig.[T] where the two 
transitions are smeared out and close together. The max- 
imum in [IqAM/AB corresponds to the narrow maxima 
seen in the field sweeps at the upper boundary of the A 
phase. It extends over a wide temperature range, since 
the transition line of the A phase for this field value dis- 
plays a very weak field dependence (cf. Fig.[l2]). We de- 



fine the transition temperatures and T^ 2 at the onset 
and the end of the maximum of hqAM/AB, respectively, 
assuming that it consists essentially of two overlapping 
maxima. As summarized below these temperatures coin- 
cide perfectly with the phase boundary inferred from the 
field sweeps. 

For lower magnetic fields the magnetization still 
displays a shallow minimum, while the maximum in 
Ho AM / AB vanishes gradually and develops into a min- 
imum for fields below 200 mT as shown in Figs.[7](b) 
through[7](d) . This corresponds still to trace (ii) in Fig. [I] 
where the two steps are now sufficiently separated. Fi- 
nally a situation more akin to trace (iv) is shown in Fig.[Tj 
where only one step is seen. We thereby continue to 
denote the lower and upper boundary of the feature in 
HoAM/AB as and T^ 2 , respectively. 

For even lower magnetic fields the magnetization de- 
velops a shallow maximum, as shown in Figs.[7](e) and 
[7](f). This corresponds to trace (iii) in Fig.[T] where the 
two steps are rather close to each other. In compari- 
son, Ho AM/ AB develops again a maximum followed by a 
minimum, where the temperature of the maximum in the 
magnetization coincides with the change of hqAM/ AB 
from a maximum to a minimum, i.e., the maximum in the 
magnetization is located in the center of the temperature 
range of the features in hqAM/AB. 

The magnetic field dependence of Re Xac does not track 
HqAM/AB in the regime of the A phase as shown in 
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FIG. 7: (Color online) Comparison of the temperature dependence of the magnetization, M, the susceptibility inferred from the 
magnetization fioAM/ AB , and the real part of the ac susceptibility, Rexac- Also shown is the ratio M/B. Data shown here 
were recorded in sample 1 for various magnetic fields in the regime of the A phase. The ac susceptibility data were recorded 
for an excitation amplitude of 0.1 mT along (100) at 10 Hz. 



MnSi, Sample 1, B||<100> 



B =0.1 mT,.-. 



CD 

CC 
CO 

< 



^ 0.8 ■ 



10Hz 
91 Hz 
911 Hz 



0.6 



< 

o 
ZL 



0.4- 




(i) 



i\w (i) n m 



J25mT 



0.12 



B =0.1 mT 

ac 



10Hz 
91 Hz 
911 Hz 



a 0.08 



- 0.04 



^ \ 215mT 
(i) (ii) = (»i) 

175mT 

" 1 1 , A 



0b 



Vx^ v ¥ tN 125mT 



26 



27 28 
T(K) 



29 



MnSi, Sample 1, B||<100> 



CD 

CC 
CO 

< 



^ 0.8 



0.6 



< 

o 

=1 



0.4 



0.12 



0O.O8- 



- 0.04 



f =10Hz 

ac 



0.1 mT 



/ A 0.3 mT 

f - 1.0 mT 



1 

(i) 

n am/ab 



(N)i(iii) 



i 

(i) 



215mT" 

(M) T..-->Wii) (iii) 

- '•* ■ W X 125mT 



f =10Hz 

ac 



(3), 



— 0.1 mT 
0.3 mT 

— 1 .0 mT 



(1) (2 ) (5)(6) m (8) 

t ^'t t ( 4 )t tt r'' 26.5 K t 
() k--- — ■- .. ~ ^ ■• • 1 



0.1 0.2 0.3 0.4 
B(T) 



FIG. 8: (Color online) Typical temperature dependence of the 
susceptibility inferred from the magnetization, /LtoAM/ AB , 
and the real and imaginary parts of the ac susceptibility, 
Re^ac and Im^ac, respectively, at selected fields for various 
excitation frequencies. Data were recorded for an excitation 
amplitude of 0.1 mT. no AM/AB was calculated from magne- 
tization data recorded for a field difference of 5 mT. Curves at 
different magnetic fields have been shifted by constant values 
for clarity. See Table |III| for a summary of the definitions of 
temperature values. 



FIG. 9: (Color online) Typical temperature dependence of the 
susceptibility inferred from the magnetization, hqAM / AB, 
and the real and imaginary parts of the ac susceptibility, 
ReXac and Im^ac, respectively, at selected fields for various 
excitation amplitudes. Data were recorded for an excitation 
frequency of 10 Hz. /j,qAM / AB was calculated from magneti- 
zation data recorded for a field difference of 5mT. Curves at 
different magnetic fields have been shifted by constant values 
for clarity. See Table |III| for a summary of the definitions of 
temperature values. 
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Fig. 7 Here we define at the point where the differ- 
ence Detween /jlqAM/AB and Rexac vanishes. Typical 
data underscoring the importance of the excitation am- 
plitude and frequency for Rexac in temperature sweeps 
are shown in Figs.|8]and[9] For increasing frequency the 
features in Re Xac become less pronounced and the peak 
in Imx ac gets slightly larger [Figs.[8](a) and[8](b)]. Like- 
wise Re Xac and Im Xac display a substantial dependence 
on the excitation amplitude, where the signal increases 
with increasing amplitude [Figs.|9](a) and [91(b)]. For all 
excitation amplitudes ReXac is smaller than fi AM/AB. 
Hence, when trying to infer phase boundaries from the 
ac susceptibility measured at a single frequency or ampli- 
tude alone there is an inherent risk for misinterpretation. 

As a final issue of the temperature dependence we 
have also explored the role of the field and tempera- 
ture history. For this purpose we determined the dif- 
ference between data recorded when heating the sample 
in an applied field after zero-field cooling (zfc-fh) and 
data recorded when cooling the sample after the field 
was applied at a temperature substantially higher than 
the helimagnetic transition temperature T c (fc). Here 
we find no differences whatsoever, as illustrated by typ- 
ical data of Re x ac and Im^ac shown in Fig.[lOj In con- 
trast, all doped systems we have studied so far, notably 
Mni-^Fe^Si, Mni-^Co^Si, and Fei-^Co^Si, display dis- 
tinct differences between zfc-fh and fc (examples will be 
given below). In fact, an observation of differences be- 
tween zfc-fh and fc in nominally stoichiometric systems 
with a complex metallurgy, say FeGe or MnGe, may hint 
at limited sample quality. 

In the spirit of a first summary we note that even 
though the temperature dependencies of M, (i AM/ AB, 
Re Xao and IniXac reported above at first sight might 
suggest a complex magnetic phase diagram in the regime 
of the A phase, the consistency of the temperature de- 
pendence with the magnetic field dependence unambigu- 
ously reveals that the phase diagram consists of a sin- 
gle pocket of the A phase. As an intrinsic property of 
the A phase the magnetization forms a tilted plateau 
with a slightly broadened first-order boundary. In fact, 
our temperature-dependent data are perfectly consistent 
with those reported by Kadowaki et gl. clearly estab- 
lishing that the phase diagram proposed in Ref.^ is in- 
correct. The complexities in the temperature dependent 
data originate in the first-order boundary of the A phase 
where the A phase coexists in a small parameter range 
with the conical state. 
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FIG. 10: (Color online) Comparison of the real and imaginary 
parts of the ac susceptibility, ReXac and ImXac, respectively, 
as recorded while field heating after zero-field cooling (zfc-fh) 
and while field cooling (fc). Data were recorded in sample 1 
for an excitation amplitude of 0.1 mT along (100) at 10 Hz. 
Sets of curves for each field value are shifted by constants for 
clarity. No differences are observed in MnSi between zfc-fh 
and fc. This is contrasted by the properties of doped systems 
such as Mni-^Fe^Si, Mni-^Co^Si, and Fei-^Co^Si, which 
display distinct differences as shown below. 



Second, is the inconsistency between the magnetization 
and ac susceptibility at the helical to conical transition 
and at the border of the A phase extrinsic or intrinsic? 

Motivated by these questions we now address changes 
of the properties reported above for samples with differ- 
ent shape and hence demagnetizing fields. We further ad- 
dress differences of the magnetization and the ac suscep- 
tibility, and thus the magnetic phase diagram, between 
magnetic field parallel to (100), (110), and (111). As a 
final aspect, we consider in further detail the frequency 
dependence of the ac susceptibility. 



A. Sample shape 



IV. FURTHER ISSUES 

The magnetization and ac susceptibility data reported 
in Sees. |III A| and |IIIB| clearly establish a single pocket 
of the A phase bounded by a regime of phase coexistence 
for field parallel to a crystallographic (100) axis. This 
raises two questions. First, what does the phase diagram 
look like for field parallel to directions other than (100)? 



Shown in Fig. 11 (a) is the real part of the ac suscepti- 
bility Re Xac as a function of applied magnetic field along 
the (100) axis for three different sample shapes. Data 
were recorded for an excitation amplitude of 1 mT and an 
excitation frequency of 10 Hz. The sample shapes studied 
here were a bar, a cube, and a platelet denoted as sam- 
ples 1, 2, and 4, respectively (cf. Table]!]). The samples 



12 



MnSi, f =10Hz, B =1 mT 



MnSi, B||<100> 



CD 

CC 







CD 

CC 




0.3 

2 

1 











27.5 K 




B II <1 00> 



Sample 1 - 6x1 x1 mm 
Sample 2-1x1x1 mm 
Sample 4 - 0.115x3x1 mm 3 




1x1x1 mm" 

-^Sample 3 - B||<111> 
-^Sample 2 - B||<110> 
Sample 2 - B||<100> 



0.1 0.2 0.3 0.4 0.5 
B(T) 

FIG. 11: (Color online) Typical ac susceptibility data for dif- 
ferent sample shapes and orientations as a function of applied 
field. Data were recorded at an excitation amplitude of 1 mT 
and 10 Hz. (a) Real part of the ac susceptibility, Rexac, for 
three different sample geometries. The pronounced qualita- 
tive and quantitative differences for different sample shapes 
are purely due to demagnetizing fields, (b) Real part of the 
ac susceptibility, Re^ac for a cubic sample shape and vari- 
ous crystallographic directions. Differences are purely due to 
differences of crystallographic orientation. 



shapes were chosen such that they differ substantially in 
terms of their demagnetizing factors. With increasing 
demagnetizing fields the characteristic transition fields 
are not only shifted to higher field values. More impor- 
tant perhaps, the signatures at the transition become less 
pronounced. Data sets recorded in sample 2 are shown 
Data as a function of field recorded 



in Figs. 14 and 15 



in sample 4, the thin platelet, are shown in Fig. 13 (the 
temperature dependent data are not shown). 

Following accurately the definitions of the characteris- 
tic field and temperature values given above, the phase 
diagrams of samples 1, 2, and 4 for field along (100) have 
been determined as a function of applied magnetic fields 
as shown in Fig. |12| As a function of temperature the he- 
limagnetic transition temperature as well as the overall 
extent of the A phase are unchanged. Not surprisingly, 
the magnetic field dependence with increasing demagne- 
tization effects is shifted to higher fields. In contrast, the 
extent of the regime of inconsistency between magnetiza- 
tion and ac susceptibility at the boundary of the A phase 
increases at the expense of the phase-pure part of the A 
phase. The increase of the regime of the inconsistency is 
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FIG. 12: (Color online) Magnetic phase diagram of MnSi as a 
function applied magnetic field along (100) for different sam- 
ple shapes. Sample 1 is a bar oriented along the applied field; 
sample 2 is a cube and sample 4 a thin platelet oriented per- 
pendicular to the applied field. With increasing demagnetiza- 
tion fields the anomalies at H c \ and for the A phase broaden 
considerably. Data points between the regime marked IM 
and the paramagnetic state at higher temperatures represent 
a crossover. 



thereby stronger than expected due to the demagnetiz- 
ing fields. This suggests that the inconsistency is not an 
intrinsic characteristic. 

Shown in Fig. [13] is the temperature and field depen- 
dence of /in d7\//d_B and Rex ac at various frequencies in 
the A phase of sample 4. Even though data values are 
rather small and therefore noisy, the same general trends 
may be observed as for the other samples. It is thereby 
also important to note that the demagnetizing effects in 
sample 4 are much stronger than for any of the other 
samples. In fact, we believe that they can no longer be 
corrected effectively by approximating the sample shape 
with an ellipsoid. 



B. Crystallographic orientation 

In order to explore systematically the importance of 
the crystallographic orientation we performed measure- 
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important difference concerns the fieid and temperature 
values of the anomalies. Clearly, as a function of tem- 
perature the extent of the A phase is much more narrow 
than for (100), in perfect agreement with the literature. 

We finally performed an analogous set of measure- 
ments for the cubic sample 3, where we studied the prop- 
erties for field parallel to (110) and (111). This way we 
first confirmed that sample 3 showed accurately the same 
field and temperature dependence for (110). Taken to- 
gether this way we obtained data that can be compared 
directly. As shown in Fig. [14] a major difference concerns 
the absolute value of the helical to conical transition. For 
field along (100) the value of B cl is largest as expected 
of the magnetically hard axis, while B c i is lowest for the 
magnetically soft (111) axis. This compares with the A 
phase, where the signatures are most pronounced for the 
(100) axis, while they are weakest for the (111) axis. 

The same trends are also seen in the temperature de- 
pendencies recorded for samples 2 and 3. For the (100) 
axis the signature of the A phase is rather pronounced, 
while is it is rather shallow and the essentially the same 
for the (110) and (111) axes with (111) being slightly 
weaker. We return to the resulting phase diagrams fur- 
ther below. 



FIG. 13: (Color online) Typical field dependence of the sus- 
ceptibility inferred from the magnetization, fiodM/dB or 
Ho AM/ AB , respectively, and the real and imaginary parts 
of the ac susceptibility, Re^ac and ImXac, respectively. Data 
shown here were recorded in sample 4 (a thin platelet) to 
track the role of large demagnetizing fields. The AC suscepti- 
bility data were recorded for an excitation amplitude of 1 mT 
along (100) at various excitation frequencies. Sets of curves 
for each temperature have been shifted by a constant for clar- 
ity. Values of the transition fields and temperatures have been 
determined in panel (a). The same values are marked in panel 
(b). 



ments on two cubic samples, denoted as 2 and 3. Starting 
with sample 2 we performed first field and temperature 
sweeps for (100) as shown in Figs. 14 (a) and 14 (b) and 



15 (a) and 15 (b), respectively. A comparison with sample 
Tithe bar for which the field was applied parallel to the 
long axis), quickly establishes that both samples display 
the same qualitative features. The cube thereby exhibits 
a slight rounding of the maxima at the various magnetic 
transitions. In contrast, there are no additional features 
between sample 1 and 2 that might hint at an important 
role other than the usual changes of the internal magnetic 
field. 

We next measured sample 2 for magnetic field applied 
parallel to the (110) axis. Since the faces of the cube 
were oriented such that the sample shape was identical 
for the two field directions, this allowed us to obtain a 
direct comparison of qualitative and quantitative differ- 
ences. Typical data are shown in Figs. 14 (c) and 14 (d) 
and 15(c) and 15(d), respectively. While the field and 



temperature dependencies are similar, the perhaps most 



C. Detailed frequency dependence 

It is finally instructive to determine the frequency de- 
pendence of the ac susceptibility in further detail. Shown 
in Fig. [16] is Rex ac in a field range crossing the A phase 
for a temperature of 27.5 K. Increasing the excitation fre- 
quency from 10 Hz to 10 kHz the maxima at the transition 
between the conical state and the A phase decrease and 
vanish. In this frequency range Rexac is always smaller 
than fig dM /dB. A similar frequency dependence is also 
seen in Imx ac , which displays a finite contribution in the 
range of the transition that decreases with increasing fre- 
quency [Fig. 16 (b)]. 

Shown in Fig. [IT] are the peak magnitudes of Rexac 
at the transition to the A phase. On a logarithmic 
frequency scale the maxima at £?ai and B^2 both de- 
crease linearly and vanish around the same character- 
istic frequency /a ~ 1.7 kHz. For frequencies exceed- 
ing /a the susceptibility at the transition approaches 
the value of the conical phase. When extrapolating the 
frequency dependence towards low frequencies the peak 
height observed in /i dM/dB would hypothetically only 
be reached for an excitation frequency around 10 -5 Hz. 

In contrast to the transition between the conical phase 
and A phase we do not observe any evidence for a max- 
imum in Re Xac at the border between the helical and 
the conical state in the frequency range between 10 Hz 
and 10 kHz that we have studied. This implies that the 
characteristic frequency / c i must be below 10 Hz; i.e., 
fci *C /a- This is consistent with spin torque exper- 
iments in the skyrmion lattice phase of MnSi, which 
demonstrate a very weak pinning of the skyrmion lat- 
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FIG. 14: (Color online) Comparison of the magnetic field dependence of the susceptibility calculated from the magnetization, 
fiodM/dB, and the real and imaginary parts of the ac susceptibility, Re^ac and ImXac, respectively, at selected frequencies. 
Data shown here were recorded for magnetic field applied along (100), (110), and (111) in samples 2 and 3. Both samples are 
cubic for the field directions studied. Thus even though data are shown as a function of applied fields they may be readily 
compared. The ac susceptibility was measured for an excitation amplitude of 1 mT. Sets of curves for the same temperature 
are shifted by constants for clarity. 



tice to defects and disorder!^ It is also consistent with 
the excellent long-range nature of the skyrmion lattice 
phase seen recently in small-angle neutron scattering. 9 
The small ac excitation may therefore be much more ef- 
ficient to generate a response of domains of the A phase 
at the border to the conical phase as compared with do- 
mains at the helical to conical transition. 



D. Relevance for other B20 compounds 

In this section we consider the question to what ex- 
tent our observations are evidence of a universal mag- 
netic phase diagram of helimagnetic B20 compounds. In 
turn this concerns foremost the disagreement between the 
magnetization and the ac susceptibility that we observe 
at the helical to conical transition and at the border of 
the A phase in pure MnSi as a generic property of all heli- 
magnetic B20 transition metal compounds. In our study 
we have addressed this question in two ways, namely 
by considering the importance of disorder in MnSi when 
doping it with Fe and, second, by considering the prop- 
erties of a rather different system - the heavily doped 
semiconductor Fei-^Co^Si. As a final aspect we discuss 
the results of ac susceptibility measurements across the 
magnetic phase diagram of FeGe.^ 

To study the importance of disorder we have measured 



simultaneously the ac susceptibility and the magnetiza- 
tion of Mn^zFezSi (x — 0.04). Previous studies have 
established that substitutional doping of MnSi with Fe 
or Co causes a suppression of the helimagnetic transition 
temperature, driving the s ystem through several quan- 
tum phase transitions! 21 * 34 ! In Mn^^FezSi with x = 0.04 
the helimagnetic transition temperature is suppressed to 
T c = 15 K. The magnetization, the ac susceptibility as 
measured at 911 Hz, and the specific heat of the same 
sample studied here were previously reported in Ref.^. 
This paper reports also the magnetic phase diagram, 
which is qualitatively reminiscent of pure MnSi. 

Shown in Fig. [18] is a comparison of the susceptibility 
calculated from the magnetization, fio dM/dB, and the 
ac susceptibility for field along (100). For the purpose 
of the study reported here field sweeps at three different 
temperatures were carried out, one well below T c at 2 K, 
the other two at 13 K and 14 K at the upper and lower 
end of the A phase, respectively. The ac susceptibility 
was thereby recorded for three different frequencies. 



Data shown in Fig. 18 (a) were obtained after zero-field 



cooling. Qualitatively these data are highly reminiscent 
of pure MnSi, with the same general features at the heli- 
cal to conical transition at B c \ and at the border of the 
A phase. In particular we observe again a clear devia- 
tion between fi dM/dB, which displays distinct maxima, 
and the ac susceptibility which does not display maxima. 
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FIG. 15: (Color online) Comparison of the temperature dependence of the susceptibility calculated from the magnetization, 
fioAM/ AB , and the real and imaginary parts of the ac susceptibility, Re^ac and Im^ac, respectively, at selected frequencies. 
Data shown here were recorded for magnetic field applied along (100), (110), and (111) in samples 2 and 3. Both samples 
are cubic for the field directions studied. Thus even though data are shown as a function of applied fields they may be 
readily compared. The ac susceptibility was measured for an excitation amplitude of 1 mT. /iqAM/AB was calculated from 
magnetization data recorded for a field difference of 5mT. Sets of curves for the same magnetic field values are shifted by 
constants for clarity. 



However, in contrast to pure MnSi, we do not observe 
any difference of the ac susceptibility between 10 Hz and 
911 Hz, the lowest and highest frequency measured, over 
the entire field range studied. This suggests that the 
disorder in slightly Fe-doped Mni_ x Fe x Si shifts the fre- 
quency dependence observed in MnSi at the border of 
the A phase to frequencies well below 10 Hz. 

The second important difference between pure MnSi 
and Mn 1 _ 2 ,Fe :c Si with x = 0.04 concerns the transition 
between the helical and the conical phase. After zero- 
field cooling B c i displays a clear temperature dependence 
and increases with decreasing temperatureP^In contrast, 
B c i is essentially temperature independent in MnSi for 
temperatures well below T c . As pointed out before, this 
suggests that disorder affects the transition at B c \ con- 
sistent with a reorientation and defect-related pinning. 
Moreover, as also pointed out before, once field cooled, 
the conical phase in Mni-^Fe^Si (a; = 0.04) is stable for 
zero magnetic field and the system does not return to 
the helical state. This conjecture was so far based on ac 
susceptibility as measured at 911 Hz only. The magnetic 
field dependence of the field-cooled state confirms these 
conclusions, as shown in Fig.[l8](b). While the anomalies 
and the difference between fiQ dM/dB and Re \&c in the 
vicinity of T c are strongly reminiscent of the zfc state, 
there is no evidence for a helical to conical phase transi- 



tion at 2 K. The reappearance of the minima at low fields 
for higher temperatures may be attributed to the effects 
of thermal excitations into the helical state. 

We turn next to the properties of Fei-zCo^Si, which 
differs from MnSi and Mni-^Fe^Si (x = 0.04) in two im- 
portant ways. First the amount of substitutional dop- 
ing is much higher, and, second, Fei-^Co^Si is a heav- 
ily doped semiconductor instead of a metal. Neverthe- 
less, the susceptibility calculated from the magnetiza- 
tion, nodM/dB, and the ac susceptibility are qualita- 
tively reminiscent of Mn 1 _ 2 ,Fe :I .Si (x = 0.04) as shown in 
Fig.[l9j Again there is no frequency dependence of the 
ac susceptibility. 

Our results shed an unexpected light on a recent study 
of the B20 compound FeGe 17 in which the magnetic 
phase diagram was inferred from the ac susceptibility 
alone. Measurements of the ac susceptibility were car- 
ried out at 1 kHz in a PPMS, i.e., an apparatus identi- 
cal to that used in our study. Moreover, a single crystal 
was studied grown by vapor transport with a rather poor 
residual resistivity ratio around The shape of the 

single crystal was ill defined making corrections of de- 
magnetizing fields difficult. Keeping these constraints in 
mind a complex phase diagram was inferred from the 
ac susceptibility, suggesting several phase pockets and 
putative evidence for mesophase formation without any 
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FIG. 16: (Color online) Comparison of the susceptibility cal- 
culated from the measured magnetization, [io dM/dB, and 
the ac susceptibility as measured at a wide range of different 
frequencies. The maxima at the lower and upper boundary of 
the A phase decrease dramatically with increasing frequency 
and vanish altogether for very large frequencies. Likewise 
maxima in the imaginary part of the susceptibility at the 
phase boundaries also decrease substantially with increasing 
frequency. 



microscopic evidence whatsoever. 

However, the qualitative field and temperature de- 
pendence of the ac susceptibility of FeGe (recorded at 
1 kHz) is strongly reminiscent of that reported here for 
MnSi and related compounds. This suggests that the 
phase diagram reported in Ref. 17 is not correct. For 
a complete phase diagram simultaneous measurements 
of the magnetization and ac susceptibility are necessary. 
Given the poor residual resistivity ratio we would thereby 
expect that the ac susceptibility does not show a fre- 
quency dependence for the same reasons that there is 
no frequency dependence in Mni-^Fe^Si (x — 0.04) and 
Fei-zCo^Si (x = 0.20). Moreover, as a vapor transport 
grown sample there may be tiny variations of the com- 
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FIG. 17: (Color online) Frequency dependence of the height 
of the peak in the ac susceptibility at the lower and upper 
boundary of the A phase denoted as Bai and Ba2, respec- 
tively. Both maxima decrease exponentially and vanish for 
frequencies above ~ 1.7 kHz and reach the size of the suscep- 
tibility in the conical phase. 



Mn Fe Si, x = 0.04, B||<100> 

1-X X ' 'II 



o 0.8 

a 

£ 0.6 

DO 

5 0.4 



T3 

i°0.2 



0.8 



cc 0.6 
ccT 

5 0.4 



T3 

i°0.2 



B =0.1 mT 

ac 
ZfC 



H dM/dB 



10Hz 
- 91 Hz 
-911 Hz 



13.0K 



B =0.1 mT 

ac 

fc 



10Hz 
91 Hz 
911 Hz 




u dM/dB 



i-utS""* 



13.0K 



0.1 



0.2 0.3 
B(T) 



0.4 



FIG. 18: (Color online) Comparison of the susceptibility 
Ho dM/dB calculated from the magnetization with the ac sus- 
ceptibility for Mni-^Fe^Si (x = 0.04). (a) Data after zero- 
field cooling, (b) Data after field cooling. The same qual- 
itative differences are observed as in MnSi with additional 
differences between the zfc and fc data in the helical state. 
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position, capable of introducing small local strains that 
might appear as several phase pockets. 



V. CONCLUSIONS 

We conclude this paper with the magnetic phase di- 
agram of MnSi as a function of internal field shown in 
Fig. |20| The phase diagrams shown here are based on 
the magnetization and ac susceptibility data recorded in 
samples 1, 2, and 3. Only the outer boundary of the A 
phase determined in sample 4 was included in Fig. 



20 



since the width of the transition surrounding the A phase 
in sample 4 is much larger than seen in samples 1,2, and 
3. We attribute this additional broadening to demagne- 
tizing fields that are more complex than can be corrected 
with the approximation as a rectangular prismP2l 

The phase diagrams shown in Fig. [20] highlight sev- 
eral key properties. Namely, the extent of the A phase 
is largest for field parallel to (100) and smallest for 
(111). This agrees with a previous report.^ Moreover, 
the lower critical field B c \ is largest for field parallel to 
(100) and smallest for (111). Both the anisotropy of the 
A phase as well as B c \ are due the same leading or- 
der magnetic anisot ropy, which is fourth order in spin- 
orbit coupling. 37 38 The same fourth-order anisotropy 
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FIG. 20: (Color online) Magnetic phase diagram of MnSi for 
field along different crystallographic directions. Phase dia- 
grams are shown as a function of internal field; the effects of 
demagnetizing fields have been corrected. Data shown here 
were inferred from samples 1 through 4. Data points between 
the regime marked IM and the paramagnetic state at higher 
temperatures represent a cross-over. 



term causes the zero-field helical order to propagate along 
(111). The temperature dependence of the upper critical 
field B C 2 is the same for all three directions independent 
of the magnetic anisotropy within the accuracy of our 
measurements. 

Our study reveals also a potential source of errors 
in the determination of the magnetic phase transitions 
in helimagnctic B20 compounds. At the phase bound- 
ary between the helical and the conical phase, as well 
as the conical and the A phase, we observe rather pro- 
nounced differences between the susceptibility calculated 
from the magnetization, fi dM/dB, and the ac suscep- 
tibility, Rexac- General thermodynamic considerations 
and the imaginary part of the ac susceptibility, Imx ac , 
identify the maxima in fi$ dM/dB as the response of large 
magnetic objects at very slow time-scales, equally consis- 
tent with a reorientation or weak first-order transition. 
Since this behavior occurs between the helical and the 
conical phase and the conical and the A phase, there is 
no reason to believe that it represents evidence for spin 
textures with nontrivial topology. 
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The disagreement between fio dM/dB and Re^ac fi- 
nally highlights the risk to determine phase transition 
lines from the ac susceptibility alone. In fact, the change 
of the isothermal magnetization inferred from the ac sus- 
ceptibility between two different field values in the conical 
phase below and above the A phase would be inconsis- 
tent for a path across the A phase and a path surrounding 
the A phase. This by itself questions strongly the phase 
diagram of FeGe reported in Ref.^. Moreover, in our 
simultaneous measurements of the magnetization and ac 
susceptibility we reproduce the magnetization reported 
by Kadowaki et aZ.™ Here our data demonstrate un- 
ambiguously that the phase diagram reported in Ref.^ 
is not correct. Instead we find that the phase diagram 
of MnSi observed in the magnetization and ac suscepti- 
bility yields a single pocket of the A phase, regardless 
of crystallographic orientation, sample shape, excitation 
frequency, and excitation amplitude of the ac suscepti- 
bility. 

We finally note that none of the rather large number 
of microscopic measurements of the A phase reported to 
date hints at a more complex magnetic phase diagram 
than reported here. Moreover, the phase boundaries 
observed in microscopic studies, notably SANS, agree 
within experimental accuracy with the phase boundaries 
we infer from the magnetization and ac susceptibility. In 



turn this questions the evidence for an exotic mechanism 
at the heart of the skyrmion lattice phaseP^ In conclu- 
sion, the combination of excellent theoretical understand- 
ing of the A phase^E^ with the similarity of our data in 
MnSi, Mni_ a; Fe 2 ;Si, and Fei-^Co^Si suggests that the 
phase diagram of MnSi represents a universal property 
of all B20 transition metal compounds with helimagnetic 
order. 

Note added in proofs. We have confirmed that the same 
differences between the magnetization and ac suscepti- 
bility reported here exist also in the magnetic phase dia- 
gram of the multiferroic insulator Cu20Se03 (cf. RefP^I). 
These results will be published elsewhere. 
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